Theory of vortex excitation imaging via an NMR relaxation measurement 
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The temperature dependence of the site-dependent nuclear spin relaxation time Ti around vor- 
tices is studied in s-wave and d-wave superconductors. Reflecting low energy electronic excitations 
associated with the vortex core, temperature dependences deviate from those of the zero-field case, 
and Ti becomes faster with approaching the vortex core. In the core region, Tj -1 has a new peak 
below T c . The NMR study by the resonance field dependence may be a new method to prove the 
spatial resolved vortex core structure in various superconductors. 

PACS numbers: 74.60.Ec, 74.25. Jb, 76.60.Pc 



Much attention has been focused on vortex physics 
both of high T c cuprates and of conventional supercon- 
ductors. Among various experimental methods, the nu- 
clear magnetic resonance (NMR) experiments!!] have been 
providing vital data in distinguishing between s-wave and 
d-wave pairing symmetries via temperature (T) depen- 
dence of the nuclear spin relaxation time Z\, which re- 
flects low-lying excitations in the superconducting state. 
The power law T-j~ oc T 3 (T 5 ) behavior is taken as defini- 
tive evidence for a line (point) node in the gap structure 
of anisotropic superconductors. This conclusion comes 
from a simple counting of the density of states (DOS) at 
the Fermi level: N(ui) oc uj(uj 3 ) for a line (point) node 
in a bulk superconductor at zero field. However, actual 
NMR experiments are performed under applied fields in a 
mixed state. Then, the-contribution of the vortex core is 
included in their datantj. Usually T\ is measured by se- 
lecting the resonance frequency at a most intensive signal 
in the resonance spectrum. However, the resonance spec- 
trum reflects information ofJnternal magnetic field distri- 
bution of the vortex latticeo. By choosing the resonance 
field, we can specify the position to detect the NMR sig- 
nal. The signal at the maximum (minimum) cutoff comes 
from the vortex center (the furthest) site. The signal at 
the logarithmic singularity of the resonance field comes 
from the saddle points of the field. By studying the po- 
sition dependence of T\ around vortices through the res- 
onance frequency dependence, we can clarify the detail 
of the vortex contribution in the NMR experiments. It 
helps us in the analysis of the standardized procedure 
extracting the gap symmetry. 

Low-lying excitation spectra around a vortex are not 
fully understood both experimentally and theoretically. 
The related problems are as follows. In the s-wave su- 
perconductors, the effect of the quantised energy level 
will appear in the quasi-particle stats3~l3. In the d-wave 
case, the low energy state around the vortex core extends 
outside the core due to the node of the superconduct- 
ing gap, and [Shows the s/H-like DOS relation (H is an 
applied field)at3. We also need to estimate the quasi- 
particle transfer between vortices (such as the path of 
the transfer and its amplitude) to study the dHvAr, os- 
cillation or transport phenomena in the mixed statefl'E3. 



The excitation around the core plays a fundamental role 
in determining physical properties of superconductors. In 
high T c cuprates, the existence or non-existence of local- 
ized core excitations in d-wave pairing case is actively 
debated. Theoretical study suggested the zero-energy 
peak in the d-wave case, instead of the quantized en- 
ergy level in the s-wave casetirli3. On the contrary, the 
scanning tunneling spectroscopy (STS) experiments ex- 
ported quantized energy level with large gap in YBCOli-3, 
and surprisingly enough-.no peak within the supercon- 
ducting gap in BSCCOOj. A part of reasons of the de- 
bate is due to limited experimental methods which di- 
rectly probe the spatially resolved core structure. So far, 
the STS was only a method to detect it. A large num- 
ber of thermodynamic or transport measurements probe 
spatially averaged quantities. Here we propose a novel 
spatially resolved means, that is, vortex imaging to see 
electronic excitations associated with a vortex core by 
using NMR, and demonstrate how the T-dependence of 
Ti is site-sensitive. Through this analysis, we are able to 
produce a spatial image of the low-lying excitation spec- 
trum around a core. A similar idea of the NMR imaging 
is actually tested, experimentally in high T c materials by 
Slichter's grpuplia and also in spin-Peiels system CuGe03 
by HorvatiaiZI. 

The position dependence of the NMR signal in the s- 
wave casfi|-was theoretically studied under some approx- 
imationsallj. Here, we calculate it microscopically from 
the wave functions obtained by self-consistently solving 
the Bogoliubov-de Gennes (BdG) equation for the ex- 
tended Hubbard model in the s- and d-wave cases. The 
eigen-energy E a and the wave functions u a {vi), v a {*i) a t 
i-site are calculated by following the method of Ref. till. 
The BdG equation for the extended Hubbard model on 
the two-dimensional square lattice is given by 
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where K ld = -tij exp[i^ f*j A(r) • dr] - Sij/j,, D i;j 



StjUAi 



TjVijAij with the on-site interaction U, the 



chemical potential fi and the flux quantum <f>o. The trans- 
fer integral t%j = t and the nearest neighber(NN) interac- 



tion Vij = V for the NN site pair r^ and r ,■ , and otherwise 



t 



, , - Vij = 0. The vector potential A(r) = ^H x r in 
the symmetric gauge. The self-consistent condition for 
the pair potential is 



Aij = --J2u a (r i )v* a (r j )ta^h(E a /2T). 



(2) 



The band filling factor (n) ~ 0.9 in our calculation. 

We consider the square vortex lattice case where near- 
est neighbor vortex is located at the 45° direction from 
the a axis. This vortex lattice configuration is suggested 
for d-w&ve superconductors, or s-waverSuperconductors 
with fourfold symmetric Fermi surfaceEHl3l£j. The unit 
cell in our calculation is the square area of N% sites where 
two vortices are included. Then, H = 2</>o/(ciV r ) 2 (c is 
the atomic lattice constant). We consider the area of N% 
unit cells. By introducing the quasi-momentum of the 
magnetic Bloch state, k = (2ir/cN r Nk)(l x ,ly) {l x ,ly = 
1, • • • , N k ), we set u a (r) = w Q (r)e lkr , v a (r) = w Q (r)e lk ' r . 
We solve Eq. (Ijj) within a unit cell under the given k. 
Then, a is labeled by k and the eigen- values obtained by 
this calculation within a unit cell. 

The periodic boundary condition is given by the sym- 
metry for the translation R = rnui + n\i2 (m and n 
are integers, ui and U2 are unit vectors of the vortex 
lattice), i.e., u Q (r + R) = u a (r)e^( r ' R )/ 2 , £ Q (r + R) = 
v a (r)e-^ (r ' R ^ 2 . Here, x(r, R) = -|jA(r) • r - nmn + 
^-(H x ro) • R in the symmetric gauge when the vortex 

center is located at r + |(ui + u 2 ). The on-site s-wave 
pair potential A a (r^) 
potential is given by 



UAa. The d x 2_ y 2-wa.ve pair 



A d (r0 = V(A £ , 



Aqj 



A-^)/4 (3) 



with A ±(M = A M±e exp[i£ / r ( ' i+ri±i)/2 A(r) ■ dr}. The 

phase factored is needed to satisfy the translational rela- 
tion A d (r) = A d (r + R)e ix ( r ' R ). 

We construct the Green's functions from E a , u a (r), 
v a (r), and calculate the spin-spin correlation function 
X-f,-(r, r',i£l n fQ. Then, we obtain the nuclear spin re- 
laxation rate, 



R(r,r') 



Im X+ ,-(r,r',in n ->n + i7))/{n/T)\n_> 
-Y,u a (r)u* at (T)[u a (r')u* a> (v')+v a (r')v* at (T')] 



xnTf\E a )6(E a -E a ,) 



(4) 



with the Fermi distribution function f(E). We consider 
the case r = r' by assuming that the nuclear relaxation 
occurs locally such as in Cu-site of high T c cuprates. 
Then, r-dependent relaxation time is given by 2\(r) = 
l/i?(r, r). In Eq. ([I|), we use S(x) = 7r _1 Im(a; — irj)^ 1 to 
consider the discrete energy level of the finite size calcu- 
lation. We typically use r\ = O.Olt. To understand the 
behavior of Ti(r), we also consider the local density of 
states (LDOS) given by 



N(r,E) 



Z^ 



[\u a (r)\'f'(E a -E) + \v a (r)\'f'(E a + E)]. (5) 



It corresponds to the differential tunnel conductance of 
STS experiments. 
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FIG. 1. The LDOS N(r, E) at sites V, A, B, C, S. (a) Position of the sites V, A, B, C, S in the square vortex lattice, where 
the nearest neighbor vortex is located in the 45° direction from the a-axis. The vortex center is shown by 0. The solid lines 
show the square atomic lattice, (b) N(r, E) in the s-wave case, (c) N(r, E) in the d-wave case. The LDOS for the S-site is 
presented by the dotted line with the solid line for the B-site. The DOS at zero-field is presented by the dotted line U with the 
solid line for the V-site. 



As for the temperature dependence of A s (r) and 
A(j(r), the vortex core radips .shrinks with decreasing T 
by the Kramer-Pesch effectli3'E3. However the shrink is 
saturated at a low temperature both in the s- and d-wave 
cases. There, the structure of A s (r) and A^(r) is almost 
independent of T. This is a quantum-limit effect which 
occurs for T/T c < Aq/Ef (Ef is the Fermi energy and 
Ao the superconducting gap at zero field. In our calcu- 
lation, Aq/Ef ~ 0.25 for d-wave, ~ 0.125 for s-wave) 
0. We calculate the low temperature behavior of Ti(r) 
by using the saturated pair potential. At higher temper- 
ature, we calculate Ti(r) by using the self-consistently 
obtained pair potential at each T. 

Figure |l| (a) shows the position of the sites (V, A, 
B, C, S) where we calculate N(E,r) and Ti(r). First, 
we see the LDOS around the vortex. The s-wave case 
(U = -2.32t, V = 0) is shown in Fig. 0(b), and the 
d-wave case (U = 0, V = — 4.20t) is shown in Fig. 0(c). 
In our calculation, N r = 20 and N^ = 8. In N(E,r) at 
the vortex center (the V-site) , the gap edge at A in the 
zero-field case (dotted line U) is smeared, and low-energy 
peaks of the vortex core state appear. In the s-wave case, 
we see some peaks above the small gap Ai (~ Aq/Ef). 
It is due to the quantization of the energy level in the 
s-wave case. In the d-wave case, the core state shows 
zero-energy peak instead of the split peaks in the s-wave 
caseO. There is no small gap. The weight of the low- 
energy states is decreased with going away from the vor- 
tex center (V^A— >B^C). Far from the vortex, N(E, r) 
is reduced to the DOS of the zero-field case. But, small 
weight of the low-energy state extending from the vortex 
core remains there. It is noted that the weight of the 
low-energy state at the S-site is larger than that of the 
B-site in the s-wave case, while the S-site is farther from 
the vortex center [see lines for the S- and B-sites in Fig 
0(b)]. It is due to the vortex lattice effect. The quasi- 
particle transfer between vortices occurs along the line 
connecting NN vortices (i.e., near the S-site). 

Next, we consider the T-dependence of Ti(r) at each 
site in Fig. [l], which reflects the LDOS discussed above. 
The NMR signal at the maximum cutoff of the resonance 
spectrum as a function of applied field or probe frequency 
comes from the vortex core at the V-site. With going 
away from the center (V^A^B^C), the resonance field 
is decreased. The signal at the minimum cutoff comes 
from the C-site. The logarithmic singularity of the res- 
onance field comes from the saddle point of the field at 
the S-site. Thus it is possible to perform the site-selective 
T\ (r) measument by tuning the resonance frequency. 

The s-wave case is shown in Fig. g. We plot Ti(r) -1 
vs. T for each site in Fig. 0(a), and re-plot it as lnTi(r) 
vs. T _1 in Fig. ||(b). We also calculate the zero-field 
case in our formulation. At the zero field, T\ ~ e A °/ T . 
Then, the slope of the InTi vs. T~ l plot gives the su- 
perconducting gap Ao, as the line U in Fig. |](b). In the 
presence of vortices, T\ deviates from the relation e A °/ T 
at low T due to the low-energy excitation around the 
vortex core. This deviation was reported in the experi- 



mentscl. In our results, reflecting the small gap Ai in the 
s-wave case, T\ shows the slope Ai at low T in the InTi 
vs. T -1 plot (see the V-site in Fig. 1(b)) as seen in Fig. 
|(b). That is, I\ ~ e Al/T . With leaving the vortex cen- 
ter, since the amplitude of the low-energy bound states is 
damped, the weight of e Al ^ T gradually decreases. Then 
the crossover temperature from e A °/ T to e Al ^ T is low- 
ered. It is noted that T\ is faster at the S-site than that 
of the B-site, while the S-site is further from the vortex 
center. This non-trivial result is due to the vortex lattice 
effect noted above. We should also notice the behavior 
of the coherence peak below T c . As seen in Fig. |](a), 
with approaching the vortex center as C— >B, the coher- 
ence peak is suppressed. But in the vortex core region 
(lines V and A), a large new peak grows at intermediate 
temperatures. This is because the LDOS at the vortex 
core has peaks at low energy Ai instead of the singularity 
of DOS at A . 
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FIG. 2. Temperature dependence of Ti(r) in the s-wave 
case at the sites V, A, B, C, S assigned in Fig. IlKa). 
(a) Ti(T c )/Ti(T) is plotted as a function of T/T c . (b) 
lnTi(T)/lnTi(T c ) is plotted as a function of T c /T. Line U 
shows the zero field case. The line N is for the normal state 
at T > T c . 
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FIG. 3. Temperature dependence of Ti(r) in the d-wave 
case at the sites V, A, B, C, S. (a) Ti(T c )/Ti(T) is plotted as 
a function of T/T c . (b) a log- log plot of (a). Line U shows the 
zero field case. The line N is for the normal state at T > T c . 
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As for the <i-wave case, we plot T^r) -1 vs. T in Fig. 
J(a), and re-plot it as a log-log plot in Fig. ||(b). At zero 
field (line U), we see the power law relation T-j - ~ T 3 of 
the rf-wave case as expected. Note that this can be seen 
only below T/T c ~ 0.1 in our case. In the presence of vor- 
tices, Ti(r)" 1 deviates from the T 3 -relation, and follows 
Ti(r) -1 ~ T at low temperature. This deviation was re- 
ported in the experiments on high-T c cupratesu. The ori- 
gin of the T-linear behavior, which is attributed to resid- 
ual density of states due to impurities or defects, is the 
low-energy state around vortices in our case. With ap- 
proaching the vortex center, the T region of the T-lincar 
behavior is enlarged and it appears from higher tempera- 
tures. As seen in Fig. [j](c) of the d-wave case, the super- 
conducting gap is buried by the low-energy state around 
vortices without the small gap of the order A^/Ep. Then, 
T-j - ~ T at low temperature in the d-wave case instead 
of the relation T\ ~ e Al / T in the s-wave case. As seen in 
Fig. y, Ti(r) -1 at the vortex center (line V) is very large 
compared with the zero-field case (line U) . It reflects the 
fact that the LDOS of the low-energy state is larger than 
the DOS of the zero- field case as seen in Fig. 0(c). This 
short relaxation may be the evidence of the low-energy 
peak in the LDOS by the low-energy core state. The 
coherence peak below T c is taken as a manifestation of 
the s-wave symmetry. In the <i-wave case, the coherence 
peak is absent. But in the vortex core region, Ty has 
a peak below T c even in the <i-wave case. We should be 
careful not to mistake this peak due to the vortex core 
relaxation as the usual coherence peak in the NMR ex- 
periment when identifying the gap symmetry. 

With increasing external magnetic field, the relax- 
ation is enhanced, because the vortex contribution is in- 
creased and the amplitude of the low energy state ex- 
tending outside the vortex core becomes large both in 
the s-wave and d-wave cases, as coinciding qualitatively 
with the observation of an orgnanic superconductor k- 
(ET) 2 Cu[N(CN) 2 ]Br by Mayaffre et alB. The details of 
the field dependence belong to a future study. 

Traditionally, the vortex contribution was considered 
as the. spin diffusion to the normal region of the vortex 
coram, and T\ is treated as the spatial average. How- 
ever, we can investigate the position dependence of Ti(r) 
around vortices through the resonance field dependence. 
This is an advantage of NMR over other methods. We 
should clarify the local mechanism of the relaxation (i.e., 
whether the relaxation occurs locally, or it is averaged by 
the spin diffusion) . It is noted that in the clean limit the 
vortex core region is not a simple core filled by normal 
state electrons! There, the characteristic T-dependence 
is expected near the vortex core other than a simple T- 
linear behavior, reflecting the rich structure of the low 
energy state around the vortex core. We expect that 
the NMR imaging study just explained here will pro- 
vide vital information for the vortex core state in high-T c 
cuprates. As for the problem whether the quantization of 
the energy levels occurs or not, T\ ~ e Al ' T if the gap Ai 



(~ Aq/Ef) is present in the excitation due to the quan- 
tization. If this small gap is absent, T-J - ~ T. As for 
the problem whether the zero-energy peak exists or not 
in the core state, the relaxation at the core becomes em- 
inently faster than that of the zero-field case (or that far 
from the vortex) at low temperature, if the zero-energy 
peak exists in the LDOS as suggested in the theoreti- 
cal study. If the peak structure is absent within Ao as 
reported in the STM experiments on BSCCO, the relax- 
ation is slow even at the vortex core as in the zero-field 
superconducting case. 

We proposed the study of the low-energy excitation 
imaging around vortices via an NMR relaxation. It may 
provide valuable information for the understanding of the 
vortex physics in high-T c superconductors as well as in 
the conventional superconductors. 

We thank M. Horvatic, K. Ishida and Y. Iwamoto for 
useful information on NMR experiments. 
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